Current constraints on Cosmological Parameters from Microwave Background 

Anisotropics. 
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We compare the latest observations of Cosmic Microwave Background (CMB) Anisotropies 
with the theoretical predictions of the standard scenario of structure formation. Assuming a 
primordial power spectrum of adiabatic perturbations we found that the total energy density is 
constrained to be Q.tot = 1-03 ± 0.06 while the energy density in baryon and Cold Dark Matter 
(CDM) are Qth^ = 0.021 ± 0.003 and n^dmh'^ = 0.12 ± 0.02, (all at 68% C.L.) respectively. 
The primordial spectrum is consistent with scale invariance, (us = 0.97 ± 0.04) and the age 
of the universe is to — 14.6 ± 0.9 Gyrs. Adding informations from Large Scale Structure and 
Supernovae, we found a strong evidence for a cosmological constant flA ~ O.TO^o os and a value 
of the Hubble parameter h = 0.69±0.07. Restricting this combined analysis to flat universes, we 
put constraints on possible 'extensions' of the standard scenario. A gravity waves contribution 
to the quadrupole anisotropy is limited to be r < 0.42 (95% c.l.). A constant equation of state 
for the dark energy component is bound to be wq < —0.74 (95% c.l.). We constrain the effective 
relativistic degrees of freedom < 6.2 and the neutrino chemical potential —0.01 < < 0.18 
and |^^,t| < 2.3 (massless neutrinos). 



m 

O 
cn 
o 

o 



I. INTRODUCTION 

The last years have been an exciting period for the 
field of the CMB research. With recent CMB balloon- 
borne and ground-based experiments we are entering a 
new era of 'precision' cosmology that enables us to use 
the CMB anisotropy measurements to constrain the 
cosmological parameters and the underlying theoreti- 
cal models. With the TOCO-97/98 ( [64], [51]) and 
BOOMERanG-97 ( [45]) experiments a firm detection 
of a first peak on about degree scales has been ob- 
tained. In the framework of adiabatic Cold Dark Mat- 
ter (CDM) models, the position, amplitude and width 
of this peak provide strong supporting evidence for the 
inflationary predictions of a low curvature (flat) uni- 
verse and a scale-invariant primordial spectrum ( [19], 
[48], [62]). 

The new experimental data from BOOMERanG 
LDB ( [53]), DASI ( [30]), MAXIMA ( [41]), CBI ( 
[57]), VSA ( [60]) and, more recently, ACBAR ( [8]), 
ARCHEOPS ( [4]) and revised and improved analy- 
sis from BOOMERanG ( [10]) and VSA ( [9]) have 
provided further evidence for the presence of the first 
peak and refined the data at larger multipole (see e.g. 
[54]). The combined data suggest the presence of a 
second and third peak in the spectrum, confirming 
the model prediction of acoustic oscillations in the 
primeval plasma and sheding new light on various cos- 
mological and inflationary parameters ( [7], [65], [59]) 



In this Rapid Communication we compare the latest 
measurements of the Cosmic Microwave Background 
Anisotropies angular power spectrum with the theo- 
retical predictions of the standard CDM scenario in 
order to constrain most of its parameters. 

Similar and careful analysis have been done re- 
cently ( [43], [57], [50]), the work presented here can 
be considered as a last-minute update of most of 
the results already published but it will also differ 
for following aspects: First of all, we will include 
the new ARCHEOPS, ACBAR, BOOMERanG and 
VSAE datasets, which provide the best determination 
to the date of region in the spectrum from the scales 
sample to COBE up to the Silk damping scales. 

Second, we will also focus on possible deviations to 
the standard scenario, like gravity waves, an equation 
of state for the dark energy wq > —1 or an extra- 
background of relativistic particles. 

Our paper is then organized as follows: In section 

II we present the analysis method we used. In Section 

III we report our results. In Section IV we discuss our 
conclusions. 



II. METHOD 

As a first step, we consider a template of adiabatic, 
A-CDM models computed with CMBFAST ( [61]), 
sampling the various parameters as follows: the phys- 
ical density in cold dark matter flcdmh^ = ^cdm = 
0.01, ...0.40, in steps of 0.01; the physical density in 
baryons il.t,h^ = ujt, ^ 0.001, 0.040, in steps of 



* However, it is important to notice that datasets ap- 
peared before April 2001 do not show presence of multi- 
ple peaks (see [66]) and that analsyses of the most recent constraints on the peak amplitude and positions (see e.g. 
datasets not based on Bayesian methods can give weaker [67]) 
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FIG. 1. Confidence contours in the Qm — i^A (left) and 
VLbh^ ~n (Right) planes from the analysis described in the 
text. 



0.001, the cosmological constant I^a = 0.0, 0.95, 
in steps of 0.05 and the curvature ilfc ~ —0.5, ....0.5 
step 0.05. The value of the Hubble constant is not an 
independent parameter, since: 



^cdm + 
1 — f^A — f^fc 



(1) 



We allow for a reionization of the intergalactic 
medium by varying also the compton optical depth 
parameter in the range Tc = 0.0, ...,0.45 in steps of 
0.05 t. 

We also vary the scalar spectral index of primordial 
fluctuations in the range ns = 0.7, 1.3 in steps of 
0.02. 

We will then restrict our analysis to flat models and, 
adding external priors as described below, wc will con- 
strain possible extensions of the standard model. In 
particular, we will consider a background of gravity 
waves, parametrized as a contribution to the CMB 
anisotropy quadrupole r = Cj/Cf . We consider the 
tensor spectral index tit to be tit — — r/6.8 for 715 < 1 
and nr = for ns > 1. 

We will also consider an equation of state for 
the dark energy wq ^ —1 sampled as wq — 
— 1.0, — 0.4 in step of 0.05. Finally we will con- 
strain an extra-background of relativistic particles, 
parametrized through an effective number of relativis- 
tic neutrinos AN^'i^ sampled as ATV^-''-^ = 0.0, 10.0 
in step of 0.5. 

For the CMB data, we use the recent results 
from the BOOMERanG-98, DASI, MAXIMA-1, 
CBI, VSAE, ACBAR and ARCHEOPS experiments. 
Where possible, we use the publicly available window 



functions and offset lognormal correction prefactors Xb 
in order to compute the theoretical band power signal 
Cb as in [12]. The likelihood for a given theoretical 
model is defined by -21n/: = {C^^ ~ C%^)Mbb' (Cf, - 
C^,) where Mbb' is the Gaussian curvature of the 
likelihood matrix at the peak. 

We include the beam and calibration uncertainties 
by the marginalization methods presented in ( [11], 
see also [69]). 

In addition to the CMB data we will also consider 
the real-space power spectrum of galaxies in the 2dF 
100k galaxy redshift survey using the data and win- 
dow functions of the analysis of Tegmark et al. ( [63]). 

To compute C'^'^^ , we evaluate pi = P{ki), where 
P{k) is the theoretical matter power spectrum and 
ki are the 49 k- values of the measurements in [63]. 
Therefore we have -2lnC'^'^^ = Y.^[Pi " {Wp)iY/dPf, 
where Pi and dPi are the measurements and corre- 
sponding error bars and W is the reported 27 x 49 
window matrix. We restrict the analysis to a range 
of scales where the fluctuations are assumed to be in 
the linear regime {k < 0.2/i~^Mpc). When combin- 
ing with the CMB data, we marginalize over a bias b 
considered to be an additional free parameter. 

Furthermore, we will also incorporate constraints 
obtained from the luminosity measurements of type 
I-a supernovae (SN-Ia). The observed apparent bolo- 
metric luminosity is related to the luminosity distance, 
measured in Mpc, by = M +5logdL{z)+25. where 
M is the absolute bolometric magnitude. The lumi- 
nosity distance is sensitive to the cosmological evolu- 
tion through an integral dependence on the Hubble 
factor d; = (1 + z) J^{dz' / H{z' ,^Q,U.M,Wq) where 
Q.Q and wq are the energy density and equation of 
state of the "dark energy" component. We evaluate 
the likelihoods assuming a constant equation of state, 
such that H{z) = poE»^i(l + z)^^+^'"^\ The pre- 
dicted m,f.f f is then calculated by calibration with low- 
z supernovae observations where the Hubble relation 
di Hqcz is obeyed. We calculate the likelihood, C, 
using the relation C ~ £o exp(— x^(^Q; ^m, wq)/2) 
where Cq is an arbitrary normalisation and is eval- 
uated using the observations of ( [26]), marginalising 
over Hq. 

Finally, we will also consider a, I — a contraint on 
the Hubble parameter, h = 0.71 ± 0.07, derived from 
Hubble Space Telescope (HST) measurements ( [25]). 

In order to constrain a parameter x we marginalize 
over the values of the other parameters y. This yields 
the marginalized likelihood distribution 



C{x)=P{x\Cb) = / C{x,y)dy 



(2) 



^We point out that values of Tc > 0.20 are in disagree- 
ment with recent estimates of the redshift of reionization 
Zrs ~ 6 ± 1 (see e.g. [27]) which point towards Tc ~ 0.05. 
However, since the reionization mechanisms is still unclear, 
we prefer to consider also greater values of Tc. 



The central values and la limits are then found from 
the 16%, 50% and 84% integrals of £(§). 
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III. RESULTS: STANDARD PARAMETERS 

In Fig. 1 wc plot the likelihood contours on the 
SIm — and fife/i^ — ns planes, using only the CMB 
data. 

As wc can sec from Figure 1 (Left Panel) the data 
strongly suggest a flat universe (i.e. ^l = ^Im + = 
1). From our CMB dataset we obtain Q = 1.03 ± 0.06 
at 95% C.L.. 

The inclusion of complementary datasets in the 
analysis breaks the angular diameter distance degen- 
eracy and provides evidence for a cosmological con- 
stant at high significance. Adding the HST constraint, 
the 2dF dataset and SN-Ia gives = 0.671^;?^, 
Oa = 0.631^ ^9 and i^A = 0.70t^;g^, all at 68% c.l.. 

Combining CMB and 2dF gives h = 0.69 ± 0.07 in 
extremely good agreement with the HST result. 

In the right panel of Fig.l we plot the CMB like- 
lihood contours in the Ofe/i^ — ns plane. As we can 
see, the present CMB data is in beautiful agreement 
with both a nearly scale invariant spectrum of primor- 
dial fluctuations, as predicted by inflation, and the 
value for the baryon density ujb = 0.020 ± 0.002 (95% 
C.L.) predicted by Standard Big Bang Nucleosynthe- 
sis (see e.g. [13]) from measurements of primordial 
deuterium. For the scalar spectral index, we found: 
Hs = 0.97 ± 0.04. However, the CMB constraint is 
also in agreement in between 2 — 0" with the lower 
BBN value 0.006 < flbh'^ < 0.017 obtained from mea- 
surements oi'^He and "^Li ( [14]) at 95% C.IL.. 

An increase in the optical depth Tc after recombi- 
nation by reionization (see e.g. [29] for a review) or by 
some more exotic mechanism damps the amplitude of 
the CMB peaks. Degeneracies with other parameters 
such as ns are present (see e.g. [5]) and we cannot 
strongly bound the value of Tc- In the range of pa- 
rameters we considered we have Tc < 0.24 at 1 — ct. 

The amount of non-baryonic dark matter is also 
constrained by the CMB data with ildmh'^ = 0.12 ± 
0.02 at 68% C.L. The presence of power around the 
third peak is crucial in this sense, since it cannot be 
easily accommodated in models based on just baryonic 
matter (see e.g. [50], [28], [47] and references therein). 

Furthermore, under the assumption of flatness, we 
can derive important constraints on the age of the 
universe to. From our cmb dataset we obtain to = 
14.6 ± 0.9 GYrs consistent with the analyses of ( [23], 
[53], [40]). 

IV. RESULTS: BEYOND THE STANDARD 

MODEL. 

As discussed before, even if the present CMB ob- 
servations can be fitted with just 5 parameters it is 
interesting to extend the analysis to other parameters 
allowed by the theory. Here we will just consider a 
few of them. 

Gravity Waves. 



The metric perturbations created during inflation 
belong to two types: scalar perturbations, which cou- 
ple to the stress-energy of matter in the universe and 
form the "seeds" for structure formation and tensor 
perturbations, also known as gravitational wave per- 
turbations. Both scalar and tensor perturbations con- 
tribute to CMB anisotropy. In most of the recent 
CMB analysis the tensor modes have been neglected, 
even though a sizable background of gravity waves is 
expected in most of the inflationary scenarios. Fur- 
thermore, in the simplest models, a detection of the 
GW background can provide information on the sec- 
ond derivative of the inflaton potential and shed light 
on the physics at ^ lO^^Gef (sec e.g. [36]). 

The shape of the Cj spectrum from tensor modes is 
drastically different from the one expected from scalar 
fluctuations, affecting only large angular scales (see 
e.g. [15]). The effect of including tensor modes is sim- 
ilar to just a rescaling of the degree-scale COBE nor- 
malization and/or a removal of the corresponding data 
points from the analysis. 

This further increases the degeneracies among cos- 
mological parameters, affecting mainly the estimates 
of the baryon and cold dark matter densities and the 
scalar spectral index ns ( [49], [38], [65], [20]). 

The amplitude of the GW background is therefore 
weakly constrained by the CMB data alone, however, 
when information from galaxy clustering and SN-Ia 
are included, an upper limit on r can be obtained. 

In Fig. 2 we plot the contraints obtained in the ns—r 
plane under the assumption of flatness and including 
the 2dF and SN-Ia data. As we can see, the possibil- 
ity of a tensor component is still in agreement with 
the combined analysis of these datasets. Including a 
conservative BBN constraint uib < 0.024 further im- 
proves the bound to r < 0.42 at 95% c.l.. Similar 
bounds have been found in a previous analysis ( [55]), 
without the VSAE, ACBAR and Boomerang revised 
datasets. 

Quintessence. 

The discovery that the universe's evolution may 
be dominated by an effective cosmological constant 
[26] is one of the most remarkable cosmological find- 
ings of recent years. One candidate that could pos- 
sibly explain the observations is a dynamical scalar 
"quintessence" field. The common characteristic of 
quintessence models is that their equations of state, 
= pI Pi vary with time while a cosmological 
constant remains fixed at wq^a = —1. Observa- 
tionally distinguishing a time variation in the equa- 
tion of state or finding wq different from —1 will 
therefore be a success for the quintessential scenario. 
Quintessence can also affect the CMB by acting as 
an additional energy component with a characteris- 
tic viscosity. However any early-universe imprint of 
quintessence is strongly constrained by Big Bang Nu- 
cleosynthesis with flQ{MeV) < 0.045 at 2a for tem- 
peratures near T ^ IMev ( [2]). 

In Figure 3 we plot the likelihood contours in the 
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FIG. 2. The 68% and 95% CL contours for the gravity 
waves contribution in the 725 — ^ — l^'i plane. 
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FIG. 3. The likelihood contours in the (JIm, wq) plane, 
with the remaining parameters taking their best-fitting 
values for the joint CMB+SN-Ia+2dF+HST analysis de- 
scribed in the text. The contours correspond to the 68%, 
95% and 99% confidence levels respectively. 



(ili/, wq) plane from our joint analyses of CMB+SN- 
Ia+HST+2dF together with the contours from the 
SN-Ia dataset only. The new CMB results improve 
the constraints from previous and similar analysis (see 
e.g., [58]), [3], [52]) with wq < -0.87 at 68% c.L. 
The current constraints are then perfectly in agree- 
ment with the Wq = — 1 cosmological constant case 
and gives no support to a quintessential field scenario 
with WQ > — 1. 

In our analysis we only consider the case of a 
constant- with- redshift wq . The assumption of a con- 
stant WQ is based on several considerations: first of 
all, since both the luminosities and angular distances 
(that are the fundamental cosmological observables) 
depend on wq through multiple integrals, they are 
not particularly sensitive to variations of wq with red- 
shift. Therefore, with current data, no strong con- 
straints can be placed on the redshift-dependence of 
Wq. Second, for most of the dynamical models on the 
market, the assumption of a piecewise-constant equa- 
tion of state is a good approximation for an unbiased 
determination of the effective equation of state 



J WQ{a)flQ{a)da 
J flQ{a)da 



(3) 



predicted by the model. Hence, if the present data is 
compatible with a constant wq = —1, it may not be 
possible to discriminate between a cosmological con- 
stant and a dynamical dark energy model. 

However one should be be very careful about draw- 
ing definitive conclusions about dark energy, since a 
constant equation of state is still an approximation of 
a real model of dark energy (see e.g. [70]). The anal- 
ysis presented here should be therefore regarded as 
a 'test' for deviations from the cosmological constant 
scenario. 

Big Bang Nucleosynthesis and Neutrinos. 

As we saw in the previous section, the SBBN 95% 
CL region, corresponding to fl^h'^ = 0.020±0.002 
(95% c.l.) (High BBN) and 0.006 < flth^ < 0.017 
(Low BBN), have a large overlap with the analogous 
CMBR contour. This fact, if it will be confirmed by 
future experiments on CMB anisotropics, can be seen 
as one of the greatest success, up to now, of the stan- 
dard hot big bang model. 

SBBN is well known to provide strong bounds on 
the number of relativistic species N^,. On the other 
hand. Degenerate BBN (DBBN), first analyzed in Ref. 
[18,24,37], gives very weak constraint on the effective 
number of massless neutrinos, since an increase in N^, 
can be compensated by a change in both the chemi- 
cal potential of the electron neutrino, = £,eT, and 
i7f,/i^. Practically, SBBN relies on the theoretical as- 
sumption that background neutrinos have negligible 
chemical potential, just like their charged lepton part- 
ners. Even though this hypothesis is perfectly justified 
by Occam razor, models have been proposed in the 
literature [1,16,17,46], where large neutrino chemical 
potentials can be generated. 
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FIG. 4. The 95% CL contours for degenerate BBN, the 
CMB results with just the age prior, combined with SN-Ia 
and combined with BBN degenerate. 



Combining the DBBN scenario with the bound on 
baryonic and radiation densities aUowed by CMBR 
data, it is possible to obtain strong constraints on the 
parameters of the model. Such an analysis was, for 
example, performed in ( [22], [42], [32], [56]) using the 
first data release of BOOMERanG and MAXIMA ( 
[6], [31]). 

We recall that the neutrino chemical potentials con- 
tribute to the total neutrino effective degrees of free- 
dom Ni, as 



= 3 H- 









(1) 


T 









(4) 



Notice that in order to get a bound on we have here 
assumed that all relativistic degrees of freedom, other 
than photons, are given by three (possibly) degenerate 
active neutrinos. 

Figure 4 summarizes the main results with the new 
CMB data for the DBBN scenario (see caption). We 
plot the 95% CL contours allowed by DBBN together 
with the analogous 95% CL region coming from the 
CMB data analysis, with only weak age prior, > 
llgyr and the 95% CL region of the joint product 
distribution C = Cdbbn'^cmb- 

We obtain the bound < 8.4, at 95% CL, which 
translates into the bounds —0.01 < < 0.24, sen- 
sibly more stringent than what can be found from 
DBBN alone. Combining CMBR and DBBN data 
with the Supernova la data [26] strongly reduces the 
degeneracy between n„i and n\. At 95% CL. we find 
AN^ < 6.2, corresponding to -0.01 < < 0.18 and 
< 2.3. 

It is however important to note that possible extra 
relativistic degrees of freedom, like light sterile neu- 
trinos, would contribute to N^, as well, and in this 
respect BBN cannot distinguish between their contri- 
bution to the total universe expansion rate and the 
one due to neutrino degeneracy. Therefore, in a more 
general framework, our estimates for Ni, can only rep- 



\Ur\ 



resent an upper bound for the total neutrino chemical 
potentials. 

Similar results have been obtained in [35], [39] and 
[33]. 



V. CONCLUSIONS 

The recent CMB data represent a beautiful success 
for the standard cosmological model. Furthermore, 
when constraints on cosmological parameters are de- 
rived under the assumption of adiabatic primordial 
perturbations their values are in agreement with the 
predictions of the theory and/or with independent ob- 
servations. 

As we saw in the previous section modifications as 
gravity waves, quintessence or extra background of 
relativistic particles are still compatible with current 
CMB observations, but are not necessary and can be 
reasonably constrained when complementary datasets 
are included. 

Since the infiationary scenario is in agreement with 
the data and all the most relevant parameters are 
starting to be constrained within a few percent ac- 
curacy, the CMB is becoming a wonderful laboratory 
for investigating the possibilities of new physics. With 
the promise of large data sets from Map, Planck and 
SNAP satellites and from the SLOAN digital sky sur- 
vey, opportunities may be open, for example, to con- 
strain dark energy models, variations in fundamental 
constants and neutrino physics. 
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